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LASER LIGHT SCATTERING DIAGNOSTIC FOR MEASUREMENT OF FLOW 
VELOCITY IN VICINITY OF PROPAGATING SHOCK WAVES 


Richard G. Seasholtz and Alvin E. Buggele 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

A laser light scattering diagnostic for measurement 
of dynamic flow velocity at a point is described. The 
instrument is being developed for use in the study of 
propagating shock waves and detonation waves in pulse 
detonation engines under development at the NASA 
Glenn Research Center (GRC). The approach uses a 
Fabry-Perot interferometer to measure the Doppler shift 
of laser light scattered from small (submicron) particles 
in the flow. The high-speed detection system required 
to resolve the transient response as a shock wave 
crosses the probe volume uses fast response 
photodetectors and a PC based data acquisition system. 
Preliminary results of measurements made in the GRC 
Mach 4, 10x25 cm supersonic wind tunnel are 
presented. Spontaneous condensation of water vapor in 
the flow is used as seed. The tunnel is supplied with 
continuous air flow at up to 45 psia; the flow is 
exhausted into the GRC laboratory-wide altitude 
exhaust system at pressures down to 0.3 psia. 

INTRODUCTION 

In this paper we describe a point measurement 
flow diagnostic being developed for high frequency 
response measurements of flow velocity. The 
instrument is being developed for use in the study of 
propagating shock waves and detonation waves in pulse 
detonation engines under development at the NASA 
Glenn Research Center. Laser based flow diagnostics 
for high frequency response measurements include Laser 
Doppler Velocimetry (LDV), which is typically based on 
scattering from particles in the one micrometer size range. 
In addition to concerns about particle lag, it is difficult to 
obtain data rates in excess of about 10 kHz. In turbulent 
flows, LDV measurements are beset by a variety of so- 
called biasing errors, caused by correlations between 
the measurement rate and flow properties. Other 
approaches, including the work presented here, are based 
on scattering from a large concentration of particles in the 
submicrometer size range. These techniques are based on 
measuring the Doppler shift of the scattered light using 
some type of frequency discriminator to convert changes 


in the optical frequency of the scattered light to changes 
in intensity. Discriminators that have been used include a 
Michelson interferometer 1 , a confocal Fabiy-Perot 
interferometer 2 and planar mirror Fabry-Perot 
interferometers 3 . Another type of optical frequency 
discriminator used for velocity measurement is based on 
atomic or molecular absorption cells 4,5 . However, this 
approach requires use of a laser that can be tuned to the 
absorption line being used. 

The approach we use is based on a previously 
developed system to measure flow density and velocity 
based on molecular Rayleigh scattering 6 . The 
frequency discriminator is a planar mirror Fabry-Perot 
interferometer used in the static, imaging mode. 
Unfortunately, molecular Rayleigh scattering using CW 
lasers in the 1-5 W range normally does not provide 
sufficiently strong scattering to allow high frequency 
response transient measurements. The alternative is to 
introduce artificial seed material into the flow to 
provide more intense scattering than that provided by 
the molecules that make up the actual flow. One 
approach for producing seed is to use a vaporization- 
condensation aerosol generator. This seed can then be 
introduced into the flow upstream of the measurement 
location. Another approach that can be used in flow 
undergoing expansion (and resulting cooling) is to rely 
on a condensable vapor. Water or carbon dioxide, for 
example, will form an aerosol by spontaneous 
condensation as the flow expands. This is the approach 
used in this paper. However, as we show, this does not 
provide sufficient scattering for high Mach number 
flows. 

Preliminary results of measurements made in the 
GRC Mach 4, 10x25 cm supersonic wind tunnel are 
presented. For these measurements, a beam from a 5 W 
CW single-frequency Nd:Vanadate laser was focused at 
a point in the flow. Scattered light was collected and 
transmitted over 130 m of optical fiber to a remote 
location, where the Fabry-Perot interferometer, 
photodetectors, and PC based data acquisition system 
were located. 
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Specific issues addressed include the design of a 
baffle assembly to allow the laser beam to pass through 
the wind tunnel test section without introducing 
excessive laser light scattering into the receiving optics. 
Also, the adequacy of using spontaneous condensation 
of water vapor as a seed material is examined. It is 
shown that spontaneous condensation of water vapor in 
ambient temperature air supplied to a relatively high 
Mach number (Mach 4) supersonic nozzle does not 
provide adequate scattering for 100 kHz sampling rate 
measurements. 


THEORY 

Light Scattering 

The frequency shift of light scattered from a moving 
object is 


fs=f-fo 


K-U 

2k 


( 1 ) 


where / 0 is the laser frequency and u is the particle 
velocity. The interaction wave vector is K = k s -ko (with 
ko and kg being the wave vectors of the incident and 
scattered light. It is convenient to introduce the velocity 
component w K = K • u/|K|, which represents the 
measured velocity component. 

The intensity of the detected scattered light can be 
expressed in terms of photoelectron count rate 
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where € is the efficiency factor (which includes the 
detector quantum efficiency), P Q is the laser power, n s is 
number density of the scattering particles, L x is the 
length of the probe volume, X is the laser wavelength, Q 
is the solid angle of the light collection optics, da/d£2 is 
the particle differential scattering cross section, % is the 
angle between the electric field vector of the (linearly 
polarized) incident light and the direction of the 
scattered light, h is Planck’s constant, and c is the 
velocity of light. For the work presented here we 
assume that sin 2 ^ = 1, which corresponds to “s” type 
scattering. 

Two factors determine the amount of scattering 
light that is available for making a measurement. One is 
the fraction of the total scattering light that can be 
collected and analyzed. This is limited by the collection 
optics and/or other optical components in the system. It 
is convenient to use the Lagrange Invariant to express 
this quantity. In our system, the optical fiber used to 
transmit the scattered light from the test cell to the 
optical processing area sets the limit on the amount of 


available scattered light. The Lagrange Invariant for the 
fiber is 


kD 

AQ = —kNA 2 (3) 

4 

where D f is the core diameter of the fiber and NA is its 
numerical aperture. For our system, D f = 0.55 mm and 
NA = 0.22, so A£2 = 0.036 mm 2 -sr. Given the f-number 
of the collection optics, we can determine the probe 
volume size using the value of the Lagrange Invariant 
set by the fiber. 


A Q = 


* 

4 4(//#) 2 


(4) 


Thus, for f/4 collection optics, the probe volume length 
L x = 0.97 mm. 

The other factor that determines the amount of 
scattered light that is available for making a 
measurement is the product of the number density of 
scattering particles n s and their differential scattering 
cross section (dcrALQ). To determine the quantity 
n s (do/dQ) we need to establish the size and 
concentration of the seed. 


Seeding 

To achieve high frequency response 
measurements, the flow must contain scattering centers 
in sufficiently high numbers so the intensity of the 
scattering will allow measurements with the desired 
uncertainty. With weak scattering, the primary source 
of error is photon statistical noise (shot noise). A 
measure of the scattering efficiency is the product of the 
number density of scattering centers n s and the 
differential scattering cross section (do/dQ). As a 
reference point, molecular Rayleigh scattering of 
532 nm light from air at NPT (normal pressure and 
temperature) has n s = 2.5xl0 25 m' 3 and (do/dQ) = 
6.1xl0' 32 mV 1 , so n s (do/dQ) = 1.5xl0" 6 mV' 1 . For 
submicron size particles, the scattering cross section is 
proportional to the sixth power of the diameter, and the 
number density of particles for a fixed amount of 
material is proportional to the inverse cube of the 
diameter. Thus the scattering efficiency is proportional 
to the cube of the particle diameter for a given amount 
of material. The optimum particle size is thus the 
largest particle that will follow the flow acceleration. 
Of course, the particles must also have a large enough 
number density so that there will be one or more 
particles in the probe volume during each sampling 
interval. 

As an example, consider the water vapor in air 
with 100% relative humidity at 20° C, which has a 
partial pressure of 0.34 psia. If we expand this humid 
air through a supersonic nozzle, most of the water vapor 
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will condense into liquid water droplets or ice particles. 
The differential scattering cross section is shown on 
figure 1 for a refractive index of 1.3 at 0.532 pm 
wavelength. Figure 2 shows the quantity n s {dol&Q) for 
Mach 2 and Mach 4 expansion. The size distribution of 
droplets formed by spontaneous condensation depends 
on the particular parameters of the expansion 7 . As an 
example, consider Mach 4 flow, and assume a particle 
diameter of 0.050 pm. Then the number density of 
scattering centers will be about n s = 0.5xl0 13 m' 3 , 
(dc/dQ) = 1.7xl0" 19 m 2 sr‘ 1 and the product n s { do/cllf ) = 
9x1 O' 7 m' 1 sr' 1 . We can then use equation 1 with P 0 = 
5 W, £ = 5%, L x = 1 mm, Q = 0.05 sr, and A = 0.532 
pm. to calculate the photoelectron count rate, N R = 3 x 
10 counts/sec. Although, this may seem to be a rather 
large rate, the number of counts in a 100 psec interval is 
only 300, and the uncertainty due the photon statistical 
noise is about 300‘ 1/2 , or 6 %. This is marginal for even 
measuring the amplitude of the signal. For measurement 
of velocity, as will be shown, we need a higher signal 
level. If the particles are smaller than the assumed 
0.050 pm, which could easily be the case, the signal 
level would lower. Note that this signal level is lower 
than the signal level for molecular Rayleigh scattering 
from air at normal temperature and pressure. However, 
for a lower Mach number flow, the number density 
would be higher. As shown in figure 2, the number 
density is about an order of magnitude larger for 
Mach 2 flow than for Mach 4 flow. 

Fabry-Perot frequency discriminator 

The velocity is determined from the scattering 
light delivered to the optical processing system shown 
in figure 3 as follows. 

Light exiting the fiber is collimated by lens L3 
(145 mm focal length) and split into two paths with an 
uncoated optical flat (BS1). About 10% of the light is 
reflected and focused by lens L4 (85 mm focal length) 
onto PMT 1 (quantum efficiency ~ 25 %). This signal 
is proportional to total scattered light. The fight 
transmitted by the beamsplitter is directed through a 
planar mirror Fabry-Perot interferometer (70 mm dia. 
mirrors, 90 % reflectivity, 11 GHz free spectral range 
(FSR), finesse ~ 15). The light exiting the 
interferometer is focused by the fringe forming lens, L6. 
This lens consists of a pair of lenses (fr2 135 mm focal 
length and f/1.2 50 mm focal length) that has an 
effective focal length of 1600 mm 

At the focal plane of the fringe pattern, the 
light is divided into two parts by a pair of mirrors 
(image dissector). A small mirror (6 mm dia.) is 
centered on the fringe pattern and directs fight from the 
central part of the fringe through lens L7 (100 mm focal 


length) to PMT 2. This small mirror is mounted on a 
larger mirror (25 mm dia.) that directs light from the 
outer part of the fringe through L8 (100 mm focal 
length) to PMT 3. The mirrors are tilted ±3° with 
respect to the optical axis. A typical calculated image 
of the inner fringe of Rayleigh scattered light is shown 
in figure 4. Note that a flow in the direction of the K 
vector, which here is in the flow direction, results in a 
positive frequency shift and increasing fringe diameter, 
while a flow in the direction opposite the K vector 
results in a negative frequency shift and decreasing 
fringe diameter. Thus, as the frequency of the Rayleigh 
scattered light increases, less fight is detected by PMT 2 
and more light is detected by PMT 3. At the mir ror 
location, the diameter of the image of the optical fiber is 
11 mm. 

Additional optics were included to provide a 
reference image of fight at the unshifted laser frequency. 
To accomplish this, several components were 
temporarily placed in the optical path using remotely 
controlled pneumatic actuators. When placed in the 
beam path, a mirror directed laser light onto a diffuser 
that scattered fight into the optical fiber. Also, a prism 
assembly (PA) was placed in the light path between the 
Fabry-Perot interferometer and the fringe fo rmin g lens. 
This served to direct light into a standard video camera. 
The video signal from this camera was digitized by a 
PC frame grabber card. A computer program analyzed 
this image and generated signals to control the Fabry- 
Perot mirror alignment. The computer data acquisition 
could be operated from two locations, one in the SWT 
control room, and the other at the remote location of the 
Fabry-Perot interferometer. 

Lower bounds for velocity measurement uncertainty 

Estimates of the measurement uncertainty for the 
technique described here, where the scattered light is 
analyzed with a planar mirror Fabry-Perot 
interferometer, are obtained by numerically calculating 
the Cramer-Rao lower bound 8 . The variance of the 
estimate of a parameter a, is given by 

v(« < )=[r- 1 ], (5) 

where no summation over repeated indices is implied. 
For Poisson statistics, T is the Fisher information matrix 
with elements 

r V 1 »K) *(»■>.) (6) 

v r(N Dq ) da, daj 
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where ) is the expected number of counts from the 

q* photodetector. The summation is over the number of 
photodetectors, and the quantities (# 0> ) depend on the 

details of the optical system and the flow parameters. 

We can calculate the Cramer-Rao lower bounds 
given by equation 5 for velocity measurements based on 
the optical configuration described in the previous 
section. This allows us to conduct parametric studies to 
determine the optimum configuration for the Fabry- 
Perot interferometer and for the light detection system. 
The expected number of photons detected in time 
interval At by detector 1 (PMT 1) is simply RN R At 
where R is the reflectivity of the beamsplitter BS1 and 
Nr is given by equation 3. The expected number of 
photons detected by PMT 2 and PMT 3 are 

<N d > = (l-R)N R Atjj S R (f)I Fp (f,e r )f L 2 dfdA (7) 

where the integrations are over frequency and the area 
of the q th detector (i.e., PMT 2 or 3), f L is the focal 
length of the fringe forming lens, and I FP is the Fabry- 
Perot instrument function given by 9 

!fp (v) - 

where iff is the phase change (neglecting any phase 
change on reflection) of the light between successive 
reflections given by 


1 + Fsin 2 [ — 
V 2 


( 8 ) 


y (/■«,) = g (9) 

c 

Here, \x is the refractive index of the medium in the 
Fabry-Perot cavity, d is the Fabry-Perot mirror spacing, 
0 r is the angle between the ray and the optic axis, and 
F = l/(sin 2 (7t/2A^) where N E is the effective finesse. In 
general, the image of a monochromatic extended source 
located in the object plane consists of a series of 
unequally spaced concentric rings. In this work, 
however, the field of view is restricted by the diameter 
of the optical fiber and includes only the inner fringe as 
shown in figure 4. 

We can readily calculate the lower bounds for 
velocity measurement uncertainty based on the detected 
light using three PMT’s as shown in figure 3. The 
uncertainties are evaluated as a function of the zero 
velocity fringe order n Q and the measured component of 
the flow velocity u K . 

It is convenient to describe the fringe location in 
terms of the fringe order rather than fringe radius because 
of the nonlinear nature of the spectral response of the 


Fabry-Perot. For example, if there is a bright fringe on 
axis, concentric bright fringes occur at integral values of 
lower orders, but the change in fringe radius decreases 
with decreasing fringe order. 

Here we define the fractional order of the fringe 
with radius r as 


For example, at the location of the image dissector, the 
fringe order of the 6 mm diameter mirror is 0.10. 
Likewise, the fringe order corresponding to the diameter 
of the image of the optical fiber is 0.33. And, if unshifted 
laser light generates a bright fringe with radius r 0 , we 
refer to this as order n 0 . Note that here we are using the 
term “fringe order” to denote change of the fringe order 
from the actual fringe order on the optical axis, which is 
2d/X. 

The change in fringe order with optical 
frequency is 

An=-^~ ( 11 ) 

FSR 

so the change of fringe order An with change of velocity 
AV is given by 


^m^m AV (12) 

FSR 

where 9 S is the scattering angle. For example, with a 
velocity of 100 m/sec, a FSR of 1 1 GHz and 120° degree 
scattering, the change in the fringe order An is 0.0296. 

Figure 5a shows the amounts of light detected from 
the inner and outer regions of the fringe pattern as a 
function of velocity for a FSR of 11 GHz, and for tie = 
0. 1 . Figure 5b shows the relation between the ratio of the 
outer/inner counts as a function of velocity. Note that this 
quantity is independent of the intensity of the scattered 
tight. Also note that this is a single valued function only 
for velocities with a magnitude less than 300 m/sec. 
The lower bound for velocity uncertainty is shown in 
figure 5c. This shows that the uncertainty is less than 
6 m/sec for velocity magnitudes less than about 
175 m/sec. It should be pointed out that the minimum 
uncertainty value can be located at any given velocity by 
shifting the value of no. Thus the important information 
given by this analysis is that the velocity measurement 
uncertainty is less than 6 m/sec over a 350 m/sec range. 

We can also examine the effect of changing the 
system parameters. For example, figure 6 shows the 
velocity uncertainty for several values of the 
interferometer free spectral range. This shows that a 
tradeoff exists between uncertainty and range. A smaller 
FSR gives lower uncertainty, but over a smaller range 
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than higher values of FSR. The results presented in this 
paper were taken with FSR =11 GHz. 


EXPERIMENT 

An experiment was conducted to evaluate the 
technique described above for measuring dynamic 
velocity. The GRC 10x25 cm supersonic wind tunnel 
was selected for these tests. The currently installed 
Mach 4 nozzle blocks were used for convenience, 
although this presented a much more demanding 
measurement situation because of the low density. 
Flows in PDE tests will generally be run at higher 
densities. 

Supersonic wind tunnel 

The GRC 10x25 cm continuous flow supersonic 
wind tunnel, configured for Mach 4 flow in the test 
section, is shown in figure 7. The air is supplied from 
the GRC laboratory-wide central compressed air system 
at pressures up to about 45 psia. The amount of water 
vapor in the air supply depended on the particular 
equipment (compressors, chillers, dryers) in use for the 
test. The flow from the tunnel exits into the laboratory 
exhaust system, which can reach pressures of less than 
0.5 psia. The mass flow for our tests was typically 
about 1 kg/sec. An electrical heater mounted in the 
tunnel air supply fine could be used to raise the 
temperature of the air supplied to the tunnel. For these 
tests, the air temperature was maintained at 25° C. A 
0.2 pm filter was located in the air supply fine to 
remove particulates. Optical access was provided by 
float plate glass windows (51 cm long, 25 cm high, 
4.8 cm thick) mounted in the tunnel. 

A set of spray nozzles was located 3 m upstream 
of the tunnel plenum to provide a capability to inject 
water aerosol to raise the humidity level and thus 
increase the quantity of scattering centers. 

An existing wedge-shaped supersonic injector 
mounted in the floor of the SWT test section was used 
to generate a moving shock wave. When the injector 
was pressurized, the resulting nozzle plume generated a 
bow shock, which moved through the probe volume as 
illustrated in figure 8. 

A 3 x 3 set of optical ports on 50.8 mm centers 
was located in the tunnel wall, as shown in figures 7 and 
9. A baffle assembly, designed to allow the laser beam 
to pass into and out of the test section without scattering 
laser light in the field of view of the receiving optics, 
could be placed in one of these ports. The 30 cm long 
assembly had a 19 mm thick window that provided the 
pressure seal. Several baffles inside the tube, as shown 
on figure 7 , were used to control the laser light scattered 
by the window. The beam then entered the tunnel 


through a final 3.2 mm diameter aperture located flush 
with the inside of the tunnel. A second similar tube was 
mounted on the opposite tunnel wall and served to 
capture the laser beam and minimize fight scattered 
back into the tunnel. 


Optical configuration 

Light from a 5W, 532 nm, single-frequency, 
Nd: Vanadate CW laser was focused by a 1000 mm 
focal length lens to a 200 jum diameter beam at the 
probe volume located in the wind tunnel test section. 
Scattered fight was collected from a downstream 
location at a 120° scattering angle and focused by a pair 
of 250 mm focal length lenses into a 0.55 mm core 
diameter, 0.22 NA, 130 m long optical fiber. The fiber 
was routed from the high-noise test cell to a quiet 
remote location, which provided a stable environment 
for the Fabry-Perot interferometer based optical 
processing system . This optical arrangement resulted 
in the measured velocity component being at 63.4° from 
the tunnel axis. The length of the probe volume was 
about 1.2 mm. Assuming that the flow is along the 
tunnel axis, the measured component = 0.44714, 
where V a is the tunnel flow velocity. 

Data acquisition 

The photoelectron pulses from the three PMT’s 
were amplified (Gain = 5) and sent to constant fraction 
discriminators (CFD). The 10 ns wide TTL level pulses 
from the CFD’s output were counted by a PC counter- 
timer board. Pulses could be simultaneously counted on 
the three channels at rates to 80 MHz indefinitely. 
Photoelectron count rates encountered in this work were 
up to 30 MHz. The accumulated counts on each channel 
were recorded at preset time intervals, typically 10 to 
1000 jusec, and stored on the computer hard disk. The 
number of counts in each time interval was given by the 
difference between adjacent values of the accumulated 
counts. 

Data processing 

The total fight scattering is measured by PMT 1. 
One approach to processing the PMT signals in order to 
obtain the velocity would be to use parameter 
estimation techniques based on the model developed for 
the calculation of the lower bounds of measurement 
uncertainty. This would entail doing a nonlinear fit of 
the data to the model function given by equation 7. 
However, because of the large quantity of data a simpler 
procedure was used. A linear relation was assumed 
between the ratio of the signals from the outer and inner 
regions (PMT’s 2 and 3). The tunnel flow velocity 
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measured using wall static pressure near the probe 
volume was used as a calibration point. The resulting 
relation between ratio of outer/inner counts to tunnel 
axial velocity is shown in figure 10. 

Results 

Data were obtained at two stations within the 
tunnel. Station 1 was 1.5 in (38 mm) downstream from 
the leading edge of the wedge-shaped injector at a 
height of 0.75 in (19 mm) above the test section floor. 
Station 2 was 3.5 in (89 mm) downstream from the 
leading edge at a height of 2.75 in (70 mm) above the 
tunnel floor. 

Figure 11a shows a time record of three channels 
of data obtained at Station 1. Data was acquired at a 
1 kHz sampling rate for 10 seconds. The axial tunnel 
velocity calculated using the linear relation between 
ratio of outer/inner signals and velocity is shown in 
figure lib. At the time the injector was pressured, the 
signal total signal level dropped from about 30 million 
counts/sec to 1.3 million counts/sec. This is probably 
due to the clean air from the injector entering the probe 
volume in place of the tunnel flow with its relatively 
large number density of scattering centers. The rms 
fluctuation in the measured velocity is about 28 m/sec 
or 4%. This can be due to either actual flow velocity 
fluctuations, shot noise, or some combination of both. 
The velocity measured in the time after the injector was 
pressurized exhibited much larger velocity fluctuations 
as a result of the much lower signal levels. 

An example of data taken at Station 2 is shown in 
figure 12. Data were acquired at a 1 kHz sampling rate 
for 30 seconds. At this location, as opposed to the data 
taken at Station 1, the level of the signal increased 
slightly as the injector was pressurized. Because the 
total signal level (7 million counts/sec without injection, 
8.5 million counts/sec with injection) did not change 
significantly, the rms fluctuation in the measured 
velocity was about the same (50 m/sec) for both states. 
This probe volume location appeared to be below and 
downstream of the bow shock, but not in the dry air of 
the injector plume. Note total count rate for this data 
was only about 25% of the count rate for the data taken 
at Station 1. This is probably due to the difference in 
the humidity level of the supply air on the two days. 

CONCLUSIONS 

The results of this initial study demonstrated a 
technique using a Fabry-Perot interferometer as an 
optical frequency discriminator to measure flow 
velocity at a point at sampling rates up to 100 kHz. 
However, obtaining acceptable measurement 
uncertainty is limited by the quantity of scattering 


centers in the flow. High frequency response 
measurements require high seeding levels. The 
measurements reported here for Mach 4 flow, where 
condensation of water vapor in the air supply provided 
the seeding, did not have sufficient strong scattering to 
give acceptable measurement uncertainties at the 
desired 100 kHz sampling rates. The introduction of 
additional moisture using water mist injection did not 
significantly increase signal levels. For lower Mach 
number expansion, the seeding levels would be 
significantly higher. For example, Mach 2 expansion 
would provide an order of magnitude higher signal 
level. Additional water vapor could also be added to 
the flow by heating the supply air and using the water 
mist injection system. 

Other flow experiments will need seed material 
provided by an aerosol generator. However, it may be 
difficult to use standard smoke generators to introduce 
seed in the high pressure flow found in many 
experiments. 

Signal levels can also be increased by use of 
higher power lasers, higher efficiency light collection 
optics and larger probe volumes. It should be noted that 
although photon counting electronics were used in this 
study, photon counting is limited to a maximum of 
about 50 MHz because of pulse pile-up errors. At 
higher signal levels, analog data acquisition should be 
used. Also, for higher signal levels, it may be 
advantageous to use high quantum efficiency 
semiconductor detectors instead of the PMT’s used 
here. 

Another approach for obtaining high sampling rate 
transient measurements is to use a pulse burst laser 
system that can provide up to a MHz pulse rate for a 
short time (up to about 100 pulses). 


REFERENCES 

Smeets, G., and A. George, “Michelson spectrometer 
for instantaneous Doppler velocity 
measurements”, J. Phys. E: Sci. Instrum, 14, pp. 
838-845,1981. 

Jackson, D.A, and D.M. Paul, “Measurement of 
supersonic velocity and turbulence by laser 
anemometry”, J. Phys. E: Sci . Instrum , 4, pp. 173- 
177,1971. 

3 Korb, C.L., B.M. Gentry and C.Y. Weng, “Edge 

technique: theory and application to the lidar 
measurement of atmospheric wind”, Appl. Opt. 31, 
pp. 4202-4213, 1992 . 

4 Crafton, J., N.M. Messersmith and J.P. Sullivan, 

“Filtered Doppler velocimetry: Development of a 
point system”, AIAA paper 98-0509, Reno 36 th , 
Jan 12-15, 1998. 


NASA/TM— 2002-21 1503 


6 



5 Kuhlman, J., P. Collins, and T. Scarberry, “Two- 

component point Doppler velocimetry data in 
circular jets," Meas. Set TechnoL , 12, pp. 395- 
408,2001. 

6 Seasholtz, R.G., J. Panda, and K.A. Elam, “Rayleigh 

scattering diagnostic for dynamic measurement of 
velocity fluctuations in high speed jets”, AIAA 39 th 
Aerospace Sciences Meeting , Reno, NV, AIAA- 
2001-0847, 2001. 


7 Lai, D.S., Generation of Monodisperse Droplets by 

Spontaneous Condensation of Flow in Nozzles , 
Ph.D. Thesis, Case Western Reserve University, 
Cleveland, OH, 1993. [Also, NASA Grant NAG3- 
831 Final Report]. 

8 Whalen, A.D., Detection of Signals in Noise , 

Academic Press, New York, pp. 324-231, 1971. 
Vaughan, J.M., The Fabry Perot Interferometer , 

History , Theory , Practice and Applications , Adam 
Hilger, Bristol, pp. 89-112, 1989. 



Fig. 1 . — Differential scattering cross section 
for scattering of 532 nm light from water 
droplets. 



Fig. 2. — Product of number density and 
differential scattering cross section for 
scattering from condensed water vapor. 




Fig. 3. — Layout of Fabry-Perot based optical processing 
system. 


Fig. 4. — Fabry-Perot interferometer fringe 
showing two regions where light is directed to 
PMT 2 (INNER) and PMT 3 (OUTER). 
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Outer/Inner count ratio 





Fig. 5. — Determination of velocity from scattering of 
532 nm laser light from condensate fog in Mach 4 
flow, (a) signals from three PMT’s; (b) ratio of 
outer/inner counts as function of measured velocity 
component; (c) lower bound for estimate of velocity 
measurement uncertainty. 



Fig. 6. — Example of lower bound for 
velocity measurement uncertainty for 
FSR = 5, 10, 15, 20, 25, and 30 GHz. 
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Supersonic Wind Tunnel With Various Configurations of Point Measurement System 
Using Supersonic Injection, Created Bow Shock Evaluation 

Station 2 - 


(6) 3/8" I.D. baffles 

(3) 1/4" I.D. baffles^' 
(3) 1/8“ I.D. baffles — 

Flipped window allows Optical access 

repositioning of point measurement system -j 

Survey probe actuator 
Config. 4(a) Mach 2 injector 
Ceiling and floor bleeds 



Optical access 3x3 array 
on 2-in. centers 

(4) 1/8" I.D. baffles on 1-in. centers 

(5) 1/4" I.D. baffles on 1-in. centers 

(2) 3/8" I.D. baffles 
on 1-in. centers 


typical\ Station 1 

Laser beers 


L — Tandem floor mounting 
for configurations 1 and 3 


Bow Shock 

1 r- Upper Umi of plum® 

! l Config. 3, 

/ / 7° half angle wedge L = .821 -in. 





Fig- 7. — Glenn Research Center 10 x 25 supersonic wind tunnel. 



Fig. 8. — Schematic view of bow shock generated by 
air injection from supersonic wedge shaped nozzle 
into supersonic cross flow. 


Fig. 9. — Photograph of baffle assembly used to reduce stray 
laser light scattering. 
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Fig. 10. — Linear calibration of ratio of 
outer/inner counts as function of SWT axial 
velocity. 





.. ... . . ' I ' , ■ , r. ■ , , ,,,,, ,, , , ,,,,,,, 

o 2 4 6 8 10 Time, sec 

Time, sec 


Fig. 11. — Data from Station 1 as injector was 
pressurized, (a) PMT signals; (b) SWT axial 
velocity. 


Fig. 12. — Data from Station 2 as injector was 
pressurized, (a) PMT signals; (b) SWT axial velocity. 
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